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Abstract 



We investigate the origin of the inter-band couphng in MgB2 by focusing 
on its unusual phononic features, namely, the strong anharmonicity of the 
phonons and the presence of both linear and quadratic electron-phonon inter- 
actions of the Su-Schrieffer-Heeger (SSH) type. The bare electronic Hamil- 
tonian has two bands with intra- and inter-band hopping, which lead to two 
decoupled hybridized bands. The phonon Hamiltonian including the anhar- 
monic terms is diagonalized approximately by a squeezing transformation, 
which causes the softening of the phonon frequency. The linear SSH coupling 
amplitude is reduced, consistently with the estimates from first-principle cal- 
culations. Additionally, the quadratic coupling generates an effective phonon- 
induced interaction between the hybridized bands, which is non-vanishing 
even in the limit of vanishing inter-bare-band hopping amplitude. 



PACS:74.20-z, 74.70-b 



I. INTRODUCTION. 



The electronic structure of the 40 K superconductor MgB2 is characterized by the pres- 
ence at the Fermi level of two hybrid bands {a and tt) of very different character [1]. This 
feature reflects itself in the experimental evidence of two gaps, which, in the absence of mag- 
netic fields, have a common critical temperature [2]. This implies an interaction between 
the (7 and tt bands contributing to the Fermi surface [3]. The microscopic origin of this 
interaction is, to the best of our knowledge, not yet clarified. The present work suggests 
that it might be traced to the unusual phononic structure of MgB2. The dominant electron- 
phonon interaction is of the Su-Schrieffer- Heeger (SSH) type, namely a modulation of the 
hopping due to the vibration of the Boron ions along the inter-site bond [4] , [5] . The unusual 



2 



features are the presence of anharmonic contributions (up to fourth order in the displace- 
ments) to the phononic Hamiltonian, and of both a hnear and a quadratic term in the SSH 
interaction [5] [6] [7]. Our aim is to give a quahtative picture of the issue by discussing a 
physically transparent model, but with no ambition of giving detailed quantitative results. 
However, we also show that our model is quantitatively consistent with the numerical results 
of Refs. [5] [6] [7] 



II. THE ELECTRONIC HAMILTONIAN. 



Our model of the electronic structure of MgB2 by a Hamiltonian has two bands, la- 
belled c and d, which hybridize through an inter-site hopping term. Then, in the real space 
representation, the bare electronic Hamiltonian reads: 



la la l(j)a l(j}a 



where J2i{j) means summing on the z nearest neigbors j of a given site /, and then on /. 
In MgB2 one expects that tfj, tff ^ tlj [1], [4]. The electron- phonon coupling parameters in 
the SSH scenario are derivatives of the hopping amplitudes with respect to the inter-site dis- 
tance. By using = N^^^'^ J2k c\a ^^P i^^^i) P^^^ to the reciprocal space representation, 
yielding : 

Hei = E (^'^ + ^tt) nl^ + E + ^tf) nt + E ^tl' {cld,^ + dlc,^) (2) 

ka ka ka 

where t^^ — z~^ Y,{j) tfj 6xp ik (jii — with x,y — c, d. To diagonalize H^i we express 

the bare operators c\^, d\^ through the hybridized operators a\^, according to: 

4a = «L COS + PL sin 4a = -«L + Pka COS (f^ (3) 

By choosing 



H^i is brought to diagonal form H^^ — Y,ka {^k'^^ta + ^k'^^ka) i '^^^^ particle energies 
in the hybridized bands given by: 
1 



+ e'^ + z {tf + 1 



(5) 



(6) 



III. THE PHONONIC HAMILTONIAN WITH ANHARMONIC TERMS. 

Following [6] , [7] we assume that the purely phononic Hamiltonian if^^ for MgB2 has to 
include, apart from the usual harmonic term, also a non-negligible quartic contribution. It 
takes the form: 

P P_ M 

where M is the Boron mass and is the frequency, at the wavevector q along the 
r — ^4 line, of the optical mode of E2g symmetry. The parameter x^p expresses the strength 
of the quartic term involving the wavevectors ±g and ±p. In MgB2, from Ref. [6], one can 
estimate x^p fa 7.8 eV"-*^. 

By quantizing the phonon field according to the usual relations: 



' (8) 



the harmonic part becomes Ml^ (^g^g + i) • When quantizing the quartic term, we 
neglect the terms with different numbers of creation and destruction operators and keep the 
remaining four-operator products only when diagonal. Namely, we approximate bLqb']jb_pb, 
{^p,q + ^p-q) ^q^-q ' whcrc 6^6^ = u^. The quartic contribution then reduces to: 



qp 



+4 E (^) (^) (1+ J + 2 1: (^) + b_A) z 



qp 



IP \ A \ A j ' ^ '"IP \ 4 / \ 4 

/ \ / QV \ /\ 



(9) 



The product u^Up is approximated in the MFA fashion, i.e. u^Up ~ ^q^^p) + {^q)^p 
(Up) (i/g) . Putting together the constant terms, we can rewrite Eq.9 as: 

E^^. {I + E-.P + (1 + ^q,-p) 



+ E^^« (^-.^S + mJE 



X 



qp 



+ const. 



Adding the harmonic contribution and defining 



X 



qp 



(10) 



(11) 



we obtain the purely phononic Hamiltonian as: 



H,^ = E ^^qXq {I + + E ^^q {^-qbl + h.,\) E ^qp + COUSt. (12) 

This form can be diagonahzed by a "squeezing " transformation [9] = 
-^qV,{bl^bl-b_^\) 



under the condition that 



tanh (2r]g) = -^E^gp ( 

q P \ 



(13) 



Notice that Eq.l3 yields 7]^ < 0. The diagonahzed Hamiltonian e^Hp^e ^ can now be 
written as: 



cosh (27?J + 2 sinh (277J E ^qp^ {KK + co^st. (14) 



By substituting rj^ from Eq.13 into Eq.l4, the renormalized frequency ur^ of the harmonic 
Hamiltonian for the squeezed phonons is written explicitly as: 



^q = ^qXq 



(15) 



' 1 - tanh" [2r]g^ 

where ^gX^ is the phonon frequency entering the quadratic part of the unsqueezed 
phononic Hamiltonian (see. Eq.l2). Thus is reduced (softened) with respect to the bare 
frequency X^fl^, in accordance with the findings of Refs. [5] [6] [7]. 
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IV. THE LINEAR ELECTRON-PHONON INTERACTION. 



The linear part of the SSH electron-phonon interaction is written, in real space, as 

where gfj — dtij/d (ui — |o = —g'ji etc., are the couphng constants. 
To obtain the Fourier-transformed form of Eq.l6 we define g^^— (1/-^) Y.ij) 9'F/6xp(iA;A^^) 
so that 



(Ok-q + €- 9-{k-g) - 9-k) = I E {sin [{k - q) ■ A^^ 

iJ) 



sm 



} = 7m (17) 



with analogous relations defining 7^^ and 7^^^. Quantization of the phonons according 
to Eq.8 leads to: 

— E l'k,q'^\a'^k-qa + lk!q^ka^k-qa + Ik^q (4:a^k-qa + ^k-qa'^ka) 

V-'V kqa 

When transformed to the hybridized fermion representation H^p' reads: 



~q I 



K) (18) 



ep 



kqa 



~ ^ E ^1 ^k%-q^ka'^k-q,a + ^k^k-qP'kcr(^k-q,(j + ^k^k-q'^k<7^k-q,a + ;,_^/3|^Q;j^._g (bLq + 

(19) 



where the effective couplings are defined as: 



^t%-q = lk,k~q COS COS + 7S-g sin (^^ sin - 27^^_g cos sin 



cd 



(20) 



= 7fe'fc-g sin <^fc sin (^^^.^ + 7^J_^ cos cos (^^.^ + 27^j^_g sin cos (^^.^^ (21) 

rS-g = Tfe^'ik-g COS if^ sin - 7^J_^ sin (/?^ cos (/p^.^ + 7^3t_<^ (cos (/^^^ cos >^^_q + sin (^^ sin ipj^_^ 

(22) 



r&-g = 7fc'ik-g sin if^^ cos (^^.^ - 7^^.^ cos (^^ sin (^^.^ - 7^^^,^ (cos (^^^^ cos (^^.^ + sin (^^ sin tf^_^^ 



(23) 



V. THE QUADRATIC ELECTRON-PHONON INTERACTION. 

According to Refs. [5] [6] [7], the electron- phonon Hamiltonian has to include also a 
quadratic term, which we write in real space as: 

Hi? = E [fiM^^^ja + fi'4A^] - ' + E /^f {die,. + c]A^) {u^ - u^) " (24) 

where fi-" = d'^tl^/d (ui — |o = /jf , etc. By defining f^'^ = fij^^^^''^ and the 

coefficients F^^^ — z {^f^'^ + f^^ — with x,y — c, d, the Fourier transform reads: 

■'■ t t dd 

E ^ka^pa'^q'^k-p-q-^kpq + ^ E ^kadpa'^q'^k-p-q-^kpq 
kpqa kpqa 



+ ^ E (dka^pa'^q'^k-p-q ^la^'ka'^qUp^k-q) ^^q (25) 
kpqa 

To pass over to the hybridized-fermion representation, let us define for short: 

Fkpq = F^pq cos cos + F^^g sin ip^ sin ip^ - F^^ sin ip^ cos ip^ (26) 

Fkpq = Kpq sin sin + F^^^ cos cos ip^ + F^^;', cos sin (p^ (27) 

^fcfq = ^kpq COS </^fe sin (^^ - F^^g sin (/^^ cos ip^ - F^^^ sin (/^^ sin ip^ (28) 

^fcpq = ^kpq sin (^jt cos (^j, - F^^g cos (^^ sin (p^ + F^^^ cos (^^ cos (p^ (29) 
Then, the quadratic coupling can be cast in the form: 

kpq 

+ ]Y E {'^pa'^kaFkpq + ^laPka^kpq + (^lal^ka^kpq + l^la^ka^kpg) U-qU-(k-p-q) (30) 
fepq 

Next, we quantize the phonons according to Eq.8. By enforcing k = p we take into 
account only the terms which do not change the number of phonons in a given mode: 



(31) 



Let us stress that our aim is to show that there are some contributions to H^^ which 
provide an effective inter-band coupUng. We do not claim to be able to treat all the terms 
in H'^: we just want to select the subset of "hot "terms. Selecting the p = q terms, then, 
Eq.30 can be written compactly as 

^ ^ E (2^.T.<. + 2<ny Ll {bl,bl + b^b_^ + + + l) 

kq 

+^ E {F^kq + F^Q («L/?fe. + l^l^ka) Ll {bl,bl + b^b_^ + i., + + l) (32) 

kq 



VI. THE ELECTRON-PHONON HAMILTONIAN IN THE SQUEEZED 

PHONON REPRESENTATION. 

Let us now introduce the squeezed phonon representation also for H'^^ + H'^ . By using 
the relation {^'^-q + ^q) = (b-q + b^^ the linear coupling term becomes: 

7^ E ^q^^'' ^k%-q'^ka'^k-q,a + ^k^k-qPkaPk-q,a + ^ k^k-q^\a(^k-q,<j + ^k,k-q(^ka^k-q,<j\ {^-q + ^ 

(33) 



Nkqa 



Then the linear coupling has a reduce amplitude, as r^^ < (see Eq.l3), consistently with 
the numerical analysis of Ref . [7] . 

For the quadratic part H^p> we get: 

e'Hi^e-' = 1 E {2F^,yta + ^Fgnl) ly^^ {bl.bl + b^b_^ + + + l) 

kq 

+ ^ E [F^kq + («L/3fe. + /3La..) Ly'^" {blqb\ + b^b_^ + + + l) (34) 

kq 
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The Eqs.33 and 34 show that the phonons induce an effective inter-band couphng 
through both the linear and the quadratic SSH interaction. However, the quadratic term 
in the second line of Eq.34 is the only one providing an inter-band coupling term which is 
non-vanishing even if {v^ — 0. It reads: 

Ly^''^ « + O («L/3.. + PUka) (35) 

where tanh {^r]^ — lim(,^^)_^o tanh (277^^ . It is interesting to discuss its behaviour in the 
limit of small hybridization of the bare bands. Prom the definition of F^^^, Ff^^^ (Eq.28 and 
29) and oi (p^{ Eq.4 ) one sees that, when the inter-band hybridization fj^ — > 0, then either 
</7j. — > 0, or — > ±7r/2. In the two cases we have, from Eqs.28 and 29: 

The reason for the minus sign in \\m.^p^-^±^^|2 F^f^ — —F^^ can be understood by recalling 
that, from Eq.3, hm^^^±^/2 4^ i^^) = ±^1^ while lim^^^±^/2 {(p^) = T^L • ^ follows: 

lim {f,% + F^Q (al(5,, + (5la,,) = -F^ (-die,, - clj,,) = (4.c,. + 4^.) 

(37) 

Hence in both cases 93^ — > 0, ±7r/2 we get the same expression for the effective inter-band 
coupling term, namely: 



= Lle'"^^ E //? {cos • A,,) - cos [{k - q) ■ A J } {die,, + ic?,,) (38) 

l{3) 

Apart from geometric factors, this term depends only on the intensity of the squeezing 
(through e^'^i) and on the amphtude of the quadratic inter-band SSH electron-phonon cou- 
pling ff^. As fij' is a second derivative of t^, it can be appreciable even if t^j itself is very 
small. Different evaluations of ff^ [1] all agree that it has an appreciable value. 
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The weakening of the hnear and quadratic electron-phonon interactions is expressed by 
the coefficients e'^" and e^''? . Their value is, in turn, set by the diagonahzation condition of 
the anharmonic phonon Hamiltonian, Eq.l3 according to: 



1 + tanh (2r]^ 



(39) 



\1- tanh (2r]g 

Therefore the squeezing effect related to the anharmonicity of the phonons also reduces 
the electron-phonon interactions. 

To conclude, let us check if the link that our model establishes between the softening 
of the harmonic frequency and the reduction of the electron-phonon coupling strength is 
consistent with the estimates of those quantities as given, e.g., in Refs. [5], [6] and [7] . The 
softening of the frequency is evaluated as 15%. [5], 25% [6] and 20% [7]. For our check 
we assume the intermediate estimate by Ref. [7]. If ujX^Q,^ — 0.80 then, from Eq.l5 , it 
follows Th (2?7g) = —0.60 and therefore e^i = 0.70, which agrees with the estimate [7] of a 
weakening of the linear electron- phonon coupling by 30%. 



VII. CONCLUSIONS. 

We have presented an analytic treatment of a two-band Hamiltonian with anharmonic 
phonons and both linear and quadratic electron-phonon interactions of the Su-Schrieffer- 
Heeger type, which should represent the essential physics of MgB2. We have shown that the 
numerical results of Ref. [7] about the phononic features of the material, namely the soften- 
ing of the effective harmonic frequency uj^ and the reduction of the linear electron-phonon 
coupling amplitude, can be interpreted as due to the phonons accomodating themselves in 
a "squeezed "state. Additionally, we have found that the quadratic electron-phonon inter- 
action generates an effective coupling between the hybridized bands due to virtual phonons. 

Acknowledgement 1 We thank J. Spalek and J. R. Iglesias for useful discussions. 
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